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Abstract 

We consider a system governed by the wave equation with index of refraction n(x), 
taken to be variable within a bounded region n C W^, and constant in \ f]. The 
solution of the time-dependent wave equation with initial data, which is localized in fi, 
spreads and decays with advancing time. This rate of decay can be measured (for d = 
1, 3, and more generally, d odd) in terms of the eigenvalues of the scattering resonance 
problem, a non-selfadjoint eigenvalue problem governing the time-harmonic solutions 
of the wave (Helmholtz) equation which are outgoing at oo. Specifically, the rate of 
energy escape from is governed by the complex scattering eigenfrequency, which is 
closest to the real axis. We study the structural design problem: Find a refractive 
index profile n*(x) within an admissible class which has a scattering frequency with 
minimal imaginary part. The admissible class is defined in terms of the compact 
support of n(x) — 1 and pointwise upper and lower (material) bounds on n(x) for 
X E Q, i.e., < n_ < n(x) < n+ < oo. We formulate this problem as a constrained 
optimization problem and prove that an optimal structure, n*(x) exists. Furthermore, 
n^(x) is piecewise constant and achieves the material bounds, i.e., n^(x) € {n_,n+}. 
In one dimension, we establish a connection between n^{x) and the well-known class 
of Bragg structures, where n{x) is constant on intervals whose length is one-quarter of 
the effective wavelength. 

1 Introduction and overview 

Many device applications, ranging from photonic to micro-mechanical require the controlled 
localization of energy within a compact region of space or "cavity". In such settings, an 
important performance-limiting loss mechanism is scattering loss, leakage from or tunneling 
out of the structure. We have in mind applications to wave phenomena in non-dissipative 
media governed by time-dependent wave equations arising, for example, in (i) electromag- 
netic waves in dielectric media, (ii) acoustic waves, and (iii) elastic waves. An important 
class of motivating examples concerns the control of light via micro- and nano-scale photonic 
crystal devices. For example, see, [JJWM08, BvFL+07]. 
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Thus, the following optimization problem naturally arises: 
Given constraints on material parameters and the size of the structure surrounding the cavity, 
how does one design a structure which maximizes the confinement time of energy? 

We next explain how the confinement-time of energy in a cavity can be expressed in terms 
of the imaginary parts of complex eigenvalue of the non- self adjoint scattering resonance 
problem (SRP). We then formulate the optimization problem, summarize the results of this 
paper, and review related work. 

1.1 Energy escape and the scattering resonance problem 

Our point of departure is the time-dependent wave equation for an inhomogeneous medium: 

n2(x) d^v{x, t) = Av{-x, t) X G R'^. (1) 

Here, n(x) denotes a spatially varying index of refraction, ^ which we assume to satisfy 
upper and lower bounds: < n_ < n(x) < < oo. We consider structures, which are 
supported in a fixed compact set, i.e., supp (n(x) — 1) = fl, where f2 is a bounded open 
subset of M^. 

Solutions to the Cauchy problem for the wave equation (1) with localized initial data 
conserve the energy: 

E[v{;t),dtv{-,t)]:= [ n2(x)|ati;(x,t)|2 + |Vt;(x,t)|Mx. (2) 

Yet, such solutions decay to zero as t — )■ oo in the pointwise or local energy sense: 

to. any compact =uteat A- CM^ / Mx.tj^x ^ 0, a. t ^ oo. 

Jk 

The rate of local energy decay or transiency of energy in a bounded set [LP89] can be derived 
by studying the solution of the initial value problem, expressed in terms of an inverse Laplace 
transform of the form 

piK+oo 

v{x,t) ~ / e^'"^^ {A + nWy^ du o vo, k > 0, 

where vq is determined by the Cauchy data at t = 0. The resolvent kernel, (A+n^ci;^)^^(x, y), 
has no poles in the upper half plane. In spatial dimensions = 1,3, it has a meromorphic 
continuation to the lower half plane, with only pole singularities. In spatial dimension 
d = 2, the resolvent kernel has a branch cut [Mel95]. These poles are called scattering 
resonances, scattering frequencies, or scattering poles. A corresponding solution is referred 

^ The wave equation in an inhomogeneous medium with position-dependent propagation speed, c(x), is 
[c~^{x)df — a) u(x, i) = 0. Let cq be the "background" homogeneous medium wave speed. The index of 
refraction, n, is defined by n(x) = co/c(x). Working in non-dimensionahzed units where cq = 1, the wave 
equation becomes (f ). 
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to as a scattering resonance mode. Other terms used are: quasi-normal mode, quasi-mode, 
or quasi-bound state. 

Due to the time-dependence e"*'^* in the inverse Laplace transform representation of the 
solution of the time-dependent initial value problem, time decay can be shown by deforming 
the contour into the lower half plane to a parallel contour along which the imaginary part 
is slightly larger than that of the scattering resonance which is closest to the real w-axis, 
i.e., the pole, w^[n], that is closest to the real w-axis gives rise to the exponential decay rate 
~ exp(- |S5a;4n]|t). See, for example, [LP89, TZOO, DC80, LVLH92, BG87] for detailed 
discussions of the role of scattering resonances and their characterizations. |53cj^[?2]| is called 
the width of the resonance and r := I^So;*!"^ is called its lifetime. 

To give a precise definition of scattering resonance solutions of the wave equation (1) as- 
sociated with the structure n(x), for which n^(x) — 1 has compact support, we first introduce 
the free space d-dimensional outgoing Green's function with pole at x: 

{ — (2ici;)~^ exp(za;|x — ?/|) d=l 
-(40-i/7«(c^|x-y|) d = 2 (3) 

(47r|x — yl)^-*- exp(za;|x — y I) = 3. 

G(|x — y\,uj) satisfies (—Ay — uj"^) Gdx — y|,u;) = (5(x — y) and, for u real, is outwardly 
radiating. 

A scattering resonance solution is a solution of the wave equation of the form u(x; ci;)e~*'^* 
which is outgoing. In particular, m(x; u) satisfies the Helmholtz equation: 

( A + cjV(x) )M(x;a;) = (4) 

Formally, writing n^(x) as n^(x) = 1 + (n^(x) — 1) and apply the outgoing Green's operator 
to (4) yields the equation: 

m(x;w)=u;M G(|x-y|,w) [^^(y) - 1] M(y;w)dy . (5) 
Jn 

A locally integrable function u(x, cu) which solves (5) is a weak solution of the Helmholtz 
equation, (4). Note that Equation (5) need only be solved for x G ^2. For x ^ fi, (5) 
gives an explicit expression for m(x) in terms of m(x) for x G ^2. We shall be particularly 
interested in bounded and piecewise constant n(x). In this case, a solution of (5) is at least 
C^(]R'^) n Hff^^iW^). The outgoing condition is encoded in u{x.]uj) being in the range of the 
outgoing Green's operator, a consequence of (5). 

Definition 1.1. 1. We refer to the integral equation (5) as the scattering resonance prob- 
lem for the structure, n(x), (SRP). 

2. The pair (w, m(x; uj)) is a scattering resonance pair if w G C and ^(x; uj) is a non-trivial 
L\^^ solution of (5). m(x, w) is called a scattering resonance mode with corresponding 
to a scattering frequency uj. 
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3. The set of scattering resonances for the structure n(x) is denoted Res„: 

Res„ := {the set of all (complex) eigenvalues, of (5) for index n(x)} (6) 

The set of scattering resonance frequencies, Res„, is discrete and lies in the open lower 
half plane, < 0. If is a scattering resonance pair, then so is {—uJ,u); the set 

ReSn is is symmetric about the imaginary axis. The set Res^ may be empty, as in the case 
where n(x) = 1 or may be non-empty, as in the explicit dimension d=l,2,3 examples in 
Appendix A. In the examples presented in dimensions 2 and 3, n^(x) is radially symmetric. 
The SRP therefore breaks into independent radially symmetric SRP's corresponding to the 
independent spherical harmonics. The scattering resonance frequency is independent of the 
particular spherical harmonic and thus we see from such examples that resonances may have 
multiplicity larger than one. 

In this paper, we study the problem of designing a refractive index profile, n{x.), subject 
to physically motivated constraints, for which there are very long-lived resonances. By the 
previous discussion, this corresponds to choosing n(x) so that there are scattering resonances 
very close to the real axis, i.e., small width, \Qu\. 

Roughly speaking, long-lived resonances can arise in the following ways: 

(A) Total internal reflection: Confinement of energy can be achieved by the mechanism 
of (nearly) total internal reflection. Consider a spherical region in 2 or 3 spatial di- 
mensions on which ?t,(x) > 1 is constant. If the "angular momentum" of the resonance 
mode is large, the mode will be strongly confined to the interface of the cavity. In 
the geometric optics approximation, the light rays have very shallow angle of incidence 
and therefore are nearly totally internally reflected. Such modes are referred to as 
whispering gallery or glancing modes and are the basis for spherical resonators; see 
section A. 

(B) Interference effects: The cavity can be surrounded by strongly reflective medium 
which is periodic of an appropriate period. In this case, wave interference effects 
provide the localizing mechanism. This is the basis for the Bragg resonator or Fabry- 
Perot cavity [JJWM08]. 

Figure 1 illustrates the difference between mechanisms (A) and (B). The left figure il- 
lustrates confinement via (nearly) total internal reflection; the refractive index is a constant 
n+ = 2 inside the circular cavity, i7 = {|x| < 1}, and n = 1 outside. Modes, /^(r)e^*^^, 
£ = 0, 1, 2, . . . (in 2D) and fi{r)Yl^{e, 0), |m| < 2£ + 1, £ = 0, 1, 2, . . . (in 3D) of increasing 
"angular momentum", i, have longer and longer lifetimes; the imaginary parts of the cor- 
responding scattering resonances tend to zero. The right figure displays radial confinement 
via interference effects; the refractive index profile consists of concentric annular regions 
alternating between = 1 and n+ = 2. 



4 





Figure 1: (left) The modulus of a mode with long lifetime {u = 4.2 — .OSSi) due to (A) total 
internal reflection because of large angular momentum {£ = 6) as in a spherical resonator, 
(right) Modulus of a mode with u = 6.5 — .039i confined by (B) interference effects, as in a 
Bragg resonator. 



Formulation of the Optimal Design Problem: We seek n(x) & A, a. specified set of 
admissible structures, having a resonance u closest to the real axis {\'^u)\ — ?■ min). Arbitrarily 
long confinement times may be achieved by either mechanism by allowing for increasingly 
large \^u\. In mechanism (A), this corresponds to rays with ever shallower angle of incidence. 
In mechanism (B), there are more wavelengths for which the wave can destructively interfere 
with itself upon multiple reflections within the cavity. We therefore impose that |3fJci;| is no 
larger than a prescribed upper bound {\^uj\ < p < oo). 

To obtain a precise formulation, we first introduce admissible sets of structures. Let 
C denote a fixed open and bounded set. Also, let <■«,_< < oo be specified. 
Then, our first admissible set is given by: 

A{fl,n_,n^) := {n: supp (n(x) — 1) C f2, n_ < n(x) < n^} (7) 

In one spatial dimension, we shall take Q = [0, L] and additionally consider the set of 
admissible structures which are symmetric: 

Asym{L,n.,n+) := {n e A{[0,L],n^,n+): n{x) = n{L - x)} . (8) 

When the choice of f2,n_,n+ are unambiguous, we simply write A and Asym- For p > 
define, for n(x) e A, Res^ C Res„: 

Res^ := {u E Res„,: \^co\ < p}. (9) 

The minimal resonance width in the set Res(^ is given by 

TP[n] := inf \^uj\. (10) 
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If Res^ = 0, then we set r^[n] = oo. We study the following 
Optimal Design Problem: 

TP(A) := inf TP\n] = inf inf (11) 

By the explicit example in section A, r^(^) < oo, for the above choices of A. 

The corresponding lifetime of the optimal resonance mode is given by t^{A) := yp\a) ' 
often used (dimensionless) quantity related to the lifetime is the quality factor (Q-factor) of 
a scattering resonance, which is defined Q := Several of the results of this paper have 

implications for the Q-factor since Q < 2rP{A) ■ 



1.2 Outline and summary of results 

In section 2, we derive variational-type identities for resonances (Proposition 2.1) in one 
spatial dimension. We use these identities to show that for structures supported on a compact 
interval, Q = [0,L], there is a general lower bound on the resonance width \'^uj\ of the 
following type: 

\^co\ > aA e-^^ (12) 

where and > depend on the constraint set, A{Q,n-,n+); see Theorem 2.2. A 
consequence of (12) is that for L sufficiently large, the maximal lifetime resonance u^, satisfies 
|3fJa;*| > 0; see Lemma 2.3. In Proposition 2.7, we also use a maximum principle argument 
to show that there exists a triangular resonance-free region in the lower half-plane. 

In section 3, we present results on the existence of a maximal lifetime resonance for 
dimensions d = 1,2, and 3. We prove, by studying the convergence of minimizing sequences, 
that there exists a structure 7T,y,(x) G A with a scattering resonance, cu^ G Res^^, satisfying 
(5) with minimal width = r^(^); see Theorem 3.1. 

In section 4 we show, in dimensions d = 1, 2, and 3, locally optimal structures Ti^(x) G A, 
are supported at the material bounds. That is, n^(x) is either or n_ for almost every 
X G fi. 

In sections 5 and 6, we specialize to the case of one- dimensional structures where we can 
prove considerably more. In section 5, we compute optimal structures using quasi-Newton 
optimization methods. These numerical observations motivate the investigations in Section 
6, of the character of 1-d optimal refractive indices, G Asym- In particular, we prove that 
optimal structures for which the associated optimal resonance is simple, are step functions, 
i.e., have a finite number of transitions between n+ and n_; see Theorem 6.4. The key to 
the proof is Proposition 6.3, a monotonicity formula for the phase of a scattering resonance 
mode. We further prove a lower bound on the interval lengths on which n = and n = n_; 
see Prop. 6.5. 

In section 6, we show that the optimal structures are related to the well-known class 
of Bragg structures, where n{x) is constant on intervals whose length is one-quarter of the 
effective wavelength. In particular, we conjecture that the lengths of the intervals, 6±, for 
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which n{x) = n± satisfies 

w r^ 1 27r 1 

d±{n^,n^, L) — > -— - — , — = -Aeff± as ra+ -> oo or L — )■ oo; 
4 \ 9iUi,\n± 4 

see Conjecture 6.3. Xes,± is called the effective wavelength. In App. C, we discuss properties 
of infinite Bragg structures, and computationally demonstrate that they optimize the spectral 
gap to midgap ratio. 

1.3 Related work 

Results on the existence of optimal scattering resonances and general bounds on the imagi- 
nary parts of scattering resonances for Schrodinger operators can be found in [Har82, HS86, 
Svi87]. Very recently, optimal designs have been considered in [Karll]. Our results for the 
Helmholtz equation make use of some of the arguments introduced in these papers. 

There has been extensive investigation during the past several decades of ideal designs for 
electromagnetic and photonic cavities, with a view toward device applications. Such design 
problems are typically formulated as an optimization problem for a particular figure of merit 
and solved using numerical optimization methods. 

The problem of maximizing the lifetime of a state trapped within a leaky cavity can be 
framed in several ways. The figure of merit can be taken to be the minimization of energy 
fiux through the boundary [LSV03] or a measure of mode localization [DS04, ABG~'~05]. In 
[KS08], the problem of minimizing \'^uj\ for a chosen resonance was investigated computa- 
tionally in both one- and two-dimensions. The 1-d problem considered here was also studied 
computationally in [HBKW08]. This work focuses on the optimization of a to minimize 
\'^uj\ which satisfies outgoing solutions of the equation dxO'dxu{x,u) + u'^n'^u^x,^) = 0. In 
particular, the variations |^ and ^ are formally computed. In [SGY06], transfer matrix 
methods were used to design low-loss 2D resonators with radial symmetry. In each of these 
papers, gradient-based optimization methods were used to solve the optimization problem. 

Genetic algorithms have also been employed to minimize energy fiux through the bound- 
ary [GPR+06, GL08]. In [EFV05, GWM02, FASKIO] the "inverse method" is employed, 
where a desired mode shape is chosen and then the material properties which produce that 
mode are found algebraically. In [BWC08], the time-dependent problem is solved to steady 
state using a finite-difference method with perfectly matched layers to approximate the out- 
going boundary conditions. The design problem is solved using a Nelder-Mead method. 

An important, related class of problems is to find photonic structures with large spec- 
tral band gaps. For the one-dimensional case, see the further discussion in Appendix C. 
Structures with optimally large band gaps have been proven to exist [CD99, Ostl2] and 
numerical methods have been applied to finding them [CDOO, BOY04, KOY05]. In [SJ03], 
topology optimization was used to find photonic crystals with optimally large bandgaps and 
also which optimally damp or guide waves. In [SH08], properties of photonic crystals with 
optimally large bandgaps are investigated. 

One property of optimal structures for (11) is that they are piecewise constant structures 
which achieve the material bounds, i.e., they are bang-bang controls. This property is also 
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realized in a number of optimization problems for eigenvalues of self-adjoint operators [Kre55, 
CM90a, CM90b, CGl+00, Ostl2] as well as for Schrddinger resonances [HS86]. In [OWll, 
Ostll] the authors consider the problem of maximizing the lifetime of a state coupled to 
radiation by an ionizing perturbation. For this class of problems, optimizers are interior 
points of the constraint set. 

Acknowledgements 

Some of the investigations of this article were motivated by discussions with P. Heider 
[HW08]. The authors wish to thank A. Barnett, D. Bindel, P. Heider, S. G. Johnson, 
R. V. Kohn, G. Ponce, J. V. Ralston, F. Santosa, O. Savin, D. Tataru, and C. W. Wong 
for stimulating interactions. Finally, we thank the referees for their helpful comments. B. 
Osting was supported in part by U.S. NSF Grant No. DMS-06-02235, EMSW21- RTG: 
Numerical Mathematics for Scientific Computing and NSF Postdoctoral Fellowship DMS- 
11-03959. M. I. Weinstein was supported in part by U.S. NSF Grants DMS-07-07850 and 
DMS-10-08855. 

2 Scattering resonances: variational identities, and 
bounds in one dimension 

In this section we derive variational-type identities and use them to obtain universal in- 
equalities for one-dimensional scattering resonances in terms of the support of n{x) — 1 and 
pointwise bounds on n(x). 

2.1 Variational identities 

In dimension one, a scattering resonance u{x, u) satisfying (5) on the domain Q = [0, L] is a 
weak solution of 



The following proposition gives a variational-type identity for the scattering resonance prob- 
lem in dimension one. 

Proposition 2.1. Let n E A, as defined in (7) be a refractive index and let {u , u{'x, u)) be 
a one dimensional scattering resonance pair on the domain Q = [0,L], i.e. a weak solution 



dlu{x;u) + co'^n^{x)u{x;uj) = 0, < a; < L 
dxu{0; oj) = —ioJu{0; u), dxu{L; u) = iuu^L; u). 



(13a) 
(13b) 
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0/ (13). Then 



^(^ ) = .L ,1 , (14a) 

Jo H-^^)\ 



/o^n2 |m(-,w)|^ 



^(^^) = — ^' \l . ' (i4b) 

Jo 

Furthermore, (14b) anc? 53w < imply: 



^ Ko..)Ph-|.(l..)P^ 

2/„''n^|«(.,u.)P 



Proof. Multiply (13a) by m(x, w) and integrate over f2 to obtain 

{^{uj^) + I'^iuj^)) [ n>|2dx= /" \u,\^d^-iu{\u{0,uj)\^ + \u{L,uj)\^) . 
Jq Jq 

Identifying real and imaginary parts yields Eq. (14). Equation (15) follows from Eq. (14b) 
and the relationship |Q=(co'^)| = 2|3f?Co'||53co'|. □ 

2.2 Lower bounds for resonances of the one-dimensional Helmholtz 
equation 

We use the variational-type identities from Sec. 2.1 to show the following universal inequality 
for one-dimensional scattering resonances. 

Theorem 2.2. Let = [0, L] C and n E A{n, n-,n+). For any scattering resonance 
u G Res„ and ^ > 0, 



\Qu)\ > min 



3exp {-{\^co 




nlL{3 + Lm^uj 





-1 , , 3exp(-4L2|SRu;|2) 

InpaMarn^>e- ^ \^^\ ^ {I + ?,nl + n\L^\^ujn ^^'^ 

Remark 2.3. In one- dimension, if p is sufficiently large, the optimal resonance, a;^, satis- 
fies I^Rcj^l > 0. This observation follows from the construction of explicit examples where 
l^*^! ^ eL(i+3n^ ) ' ^^^^ includcs simple piecewise constant structures [HW08] or numerically 
constructed refractive indices (see Section 5). 

We follow the strategy of [Har82] to prove Theorem 2.2. Theorem 2.2 relies on the 
following two lemmata, which we shall prove first. In this section only, we normalize the 
resonance state by assuming, without loss of generality, 

1 = u(0) < |u(L)|2 (17) 
(otherwise we make the substitution x ^ L — x). 
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Lemma 2.4. Let u{x; u)),uj denote a scattering resonance pair for the one- dimensional scat- 
tering resonance problem, (13), defined for < x < L. Then, we have the pointwise bound 

\u{x]uj)\ < a/1 + Iw^x^ exp ( / {x — y)n^{y)dy 



Proof. Assuming (17), the boundary condition (13b) is written m(0) = 1, ^^(O) = —iuJu{Q) 
and Ux{L) = iuu{L). Integrating twice, we obtain the integral equation 



X fy 

2/ 



u{x) = 1 — iujx — uj / / n {z)u{z)dz<ly 



^0 



Integrating the outer integral by parts, we obtain \j^j^n^{z)u{z)(iz(iy\ < J^(x—y)n'^{y)\u{y)\dy 
and thus \u{x)\ < a/1 + Iw^x^ + lw^ jQ{x — y)n'^{y)\u{y)\ dy. Equation (18) now follows from 
Gronwall's inequality. □ 

Lemma 2.5. Assuming the same hypotheses as in Lemma 2.4, 

[ n>p dx < nlL exp (Iwpn^L^) (l + \uj\'^L^/3) (19) 
Jo 

where n+ = maXa;g(o,L) n(x). 

Proof. Using Lemma 2.4, we compute 

n^|u|^dx< J n^(x)(l + |a;|^x^) exp ^2|a;|^ y (x — y)n^{y) dy^ dx 



<nlexp(^2\u\^ {L - y)n\y) dy^ j ^ 
<<exp(|a;|X^2) (L + |u;|2lV3) 



1 + IwPx^ dx 



as desired. □ 
Proof of Theorem 2.2. Using (15), (17), and Lemma 2.5, we compute 



exp (— 


\uj 




nlL{l + \ 


uj\ 


|2L73) 



exp (— ( 


? + \ 






n^L(l + (|3fJu;|2 + 


So; 


|2)L2/3) 



This is a nonlinear inequality for \'^uj\. Note that /(x) is a monotonically decreasing function 
for X > with / | as x t cso. Thus, for ^ > \'^uj\, /(^) < /(|53a;|) < \'^uj\. Thus for all 
,^ > 0, > min[,^, /(.^ )]. To obtain the optimal bound, one would choose ^ = such 
that ,^0 = /('Co)- For simplicity, we choose = {.'<^+L)~^ ■ If > e"^, we find that 
min[eo, /(Co)] = /(^o) for all ^uj and (16) follows. □ 
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Remark 2.6. Theorem 2.2 also gives an upper bound for the quahty factor, defined Q :- 



. In particular this bound shows that Q J, as JJa; J, 0. 



|SRa;| 
2|Sw| • 

The following proposition shows that there is a triangular resonance-free region in the 
lower-half complex plane. 

Proposition 2.7. Let n G Asym, i-^-, n_ < n{x) < n+ and n{L — x) = L{x). If u E ReSn 
is a one-dimensional Helmholtz resonance satisfying (13) with d = 1, then u ^ {u: \'^u)\ > 
\^u\ and \'^oj\ < :;^}- 

Proof. The proof of this theorem follows [HS86]. Let \Quj\ > \^u}\ and we'll show that 
\Quj\ > Using Eq. (15) and \u{0)\ = 1^(1^)1 (see Prop. 6.1), we have 

l^'^l > (21) 



f \u\^dx 



Kato's inequality [RS79] and \^uj\ > \^uj\ ^{lu^) < then give A\u\ > ^ (^Am 

—^{u'^)n'^\u\ > 0. We now apply the maximum principle to the subharmonic function |ti(a;)| 
to obtain \u{x)\ < |'u(0)|. It now follows from Eq. 21 that > □ 

3 Existence of a solution for the spectral optimization 
problem 

In this section, we consider the spectral optimization problem in dimension d = 1,2,3 with 
admissible set A{fl, n_, n+), as defined in Eq. (7), the set of n(x) satisfying upper and lower 
bounds on the compact set Q with n(x) = 1 for x ^ fi, i.e. 

n_ < n{x) < n^, x E Q, and n{x) = 1, x ^ Q. 

Recall Res^ as defined in (9) is the set of scattering frequencies, u, for the structure n(x) 
such that \^u\ < p. 

Theorem 3.1. Consider the scattering resonance problem on M"^, d = 1,2,3. Fix p > 0. 
Assume that there exists n E A such that Res^ ^ 0. Then the double infimum, defined in 
(11), is strictly positive and is attained for an admissible structure. That is, there exists 
G A, with associated longest-lived resonance mode, u^,, of frequency G Res^ and such 
that \^uj^\= TP{A) > 0. 

Proof. Since Res^ 7^ and Res^ C {u: < }, we have < \^uj\ < 00 and there is a 
minimizing sequence {^mj^^i C A, such that 

inf{ l^^wl : w G Res^^ } > 0, as m t 00 . (22) 

We first show that is attained and then conclude the proof by showing > 0. 
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Let ( 0Jm,Umi^,(^m) ) denote a sequence of resonance pairs corresponding to this min- 
imizing sequence of structures in A. Since [— p,p] is compact, there exists a convergent 
subsequence, which we continue to denote by {um}, with Um <^*- 

Since Q is bounded, we have fl C -B_r(0), where Br{0) denotes the open ball of radius R 
about the origin. By linearity and boundedness of -Br(O) we can impose the normalization 
||^m||L2(Bfl(o)) = 1- Squaring the differential equation for Um, —^u^^x) = ci;^n^(x)Mm(x), 
and integrating over B^{0), we obtain 

\AumMf dx. < luml'^ n\ I I Mm, fx) I ^(ix . 



Br{0) JBr{Q) 

Therefore, since ||Mm.||L2(5^(o)) = 1; the sequence is uniformly bounded in H^{0), 

for any open O, whose closure is a compact subset of Bii{0). Thus, by Rellich's Lemma, 
for any s < 2 there exists G H'^{Q) and a strongly convergent subsequence converging 
to M^. Moreover, Um is uniformly Holder continuous with exponent a G (0,1/2). Thus, 
{umi'x.)} is uniformly bounded and equicontinuous. Therefore, there exists a subsequence, 
again denoted {um}, such that Um u^, uniformly in VL. It follows that ||u*||L2(f^) = 1 and 
thus is nonzero. 

The uniform bounds n(x) C n+] imply that the sequence {rim} is uniformly bounded 



I.e., 



in L^(f2), and therefore along some subsequence converges weakly in L^(i7) to some n- 
rim n^, with G A. Furthermore, nm(x) — ?• n*(x) a.e. in Vt. 

It remains to show that (m*(-, o;^), o;^) is a resonance pair, i. e., u^, is non-trivial (established 
just above) solution of (5). For each m G N, we also have that 

Um{^)=ul^j G'(|x-y|,a;™)[n^(y) - l]nm(y)dy. (23) 
Jn 

As shown above, the right hand side of (23) converges to Mjf(x) uniformly on f2. Therefore 
to establish (5), it suffices to show that 



I G'(|x-y|,a;J[n^(y)-lK„(y)dy ul I G(|x-y|, u;,)[n2(y) - lK(y) dy (24) 



for X G For each x G ^2, the integrand of (24) converges in y pointwise a.e. to the 
expression with the subscript m replaced by -k. Also for each x G fi, the integrand (as a 
function of y) is dominated by an integrable function; there exists a constant K such that 



|G(|x-y|,w„)[n^(y) - l]M„(y)| < fsT |x-y| 



2-d 



uniformly in m. Hence, by the dominated convergence theorem, (24) holds. 

Finally, we claim that > 0. Suppose not. Then F^ = and the scattering resonance 
problem (5) has a non-trivial solution with real frequency w^. By the unique continuation 
principle [CK98, KT06] = 0, a contradiction. □ 
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Remark 3.2. In the example of Appendix A, we showed that in dimensions d = 2,3 for 
Q = {|x| < a} and n{x) = 1 + rioln, the resonances approach the real axis as the angular 
momentum £ f oo. In section 5, we show that in dimension d = 1, the sequence for 
increasing p'^ cx3 is such that \'^uj'^\ i with |3fJa;^| ^ ^ oo. Thus for d = 1,2,3 an 
optimal solution for the limit p t is not achieved. This result contrasts the behavior of 
optimal resonances of the Schrodinger operator [Har82, HS86, Svi87]. 

4 Local optimizers are piecewise constant structures 
which saturate the constraints 

In this section, we focus on properties of locally optimal solutions of the design problem 

min Vln]. (25) 

We begin by computing the variation of u with respect to changes in the index ?t,(x) G 
A C L°°. Assume a mapping J' : A C, n ^ J[n\. We say the mapping J is Frechet 
differentiable at n G ^ if there exists a mapping n G ^ t-^ G Li(fi) such that for all 
5n G L°° such that n + 5n & A and ||5n||Loo sufficiently small we have 

J[n + 5n]- J[n]- l^^,5n) 0, as ||(5n||Loo(n) ^ 0. 



Here (/, g) = /(x) g{x) dx, defined for f E and g G L°°. 

Proposition 4.1. Let (cj, ^(x,^)) be a nondegenerate scattering resonance pair of the scat- 
tering resonance problem (5) for index of refraction, n{x). Then 

1. The first variation of uj[n] with respect to n{x.) is given by 

-^(x) = -2a n(x) M^'^ (26) 
on 

where a G C \ {0} depends on u. 

2. In one- dimension, with Vt = [0,L], the first variation is given by (26) with 



' = 2u /"'' nV + z[u\0) + u\L)] = - ^ ul + u?n^v? (27) 



a 



Proof. The proof of Proposition 4.1 is given in Appendix B. □ 

Remark 4.2. The following are properties of the variation |^(x): 

1. From the explicit expression for |^(x), we see that it is invariant under 'u(x) i— )■ 
cm(x), c G M. 
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2. For finite degeneracy, (26) remains valid, except the correct (perturbation dependent) 
scattering resonance mode must be chosen from the eigenspace. 

3. Suppose d = 1 and = [0,L]. If n{x) is symmetric, i.e., n[x) = n{L — x), then j^{x) 
is also symmetric. This follows from Prop. 6.1. 

4. If(i = 2or(i = 3, and n(x) is radially symmetric, then |^(x) is not radially symmetric 
if the mode m(x, uj) has a nontrivial angular dependence. 



Proposition 4.3. Let n^{x.) E A be a local minimizer of the design problem (25). If OJi, 
is nondegenerate and \^uj^\ < p, then n^(x) is piecewise constant and attains the material 
bounds, i.e., [^^,t(x) — nJ\ [^*(x) — = for almost every x G 



Remark 4.4. Computational evidence for Proposition 4.3 has been reported on in [KS08, 
HBKW08]. A similar result in one-dimension with slightly different boundary conditions 
is given in [Karll]. This phenomena has also been studied in self-adjoint systems [Kre55, 
CM90a, CM90b] and for the analogous problem for Schrodinger resonances [HS86]. In control 
theory, n{^) would be referred to as a "bang-bang control." 

Remark 4.5. Proposition 4.3 has significance in applications because 

1. the optimization problem is reduced to finding the (ci— l)-dimensional interface between 
regions with constant n(x), and 

2. the computation of Helmholtz resonances for piecewise constant n(x) can be performed 
more efficiently using methods which utilize the free-space Green's function (3), e.g., 
meshless methods, the method of particular solutions, and layer potential methods, 

3. standard manufacturing techniques for photonic crystals involve creating air holes in 
a dielectric media [JJWM08]. Thus, in practice, ?t,(x) only takes two values. 



Proof. Let A dVthe defined to be the set A := {x G ri: n_ < ^^.(x) < n+} and B d Ahe 
an arbitrary set. Local optimality requires that 



— , Ifi 
on 



$5 (^ul j ui^ = $5 (aculuK^)) =0 Vx G A 



since B is arbitrary. We conclude that au^uKx.) is a real- valued function on A. 

In the neighborhood of any x such that Uj,(x) ^ 0, we can choose the sign such that 
t'(x) := a/ ±aujlul{x.) is a real- valued function which satisfies the complex Eq. (4): Av + 
3fJ(a;^)n^f = —iQ{ul)n'^v. Since the left-hand side of this equation is purely real and the 
right-hand side is purely imaginary, we conclude that v = 0. Thus u^, = on A. 

Since m(x) satisfies (4), we know that by the unique continuation principle [CK98, KTOl], 
it cannot vanish on an open set of Q, or else m*(x) = 0, a contradiction. Thus meas(A) = 
and we identify it with the zero level set of ^3 (auluKx.)) . □ 
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Proposition 4.3 is further strengthened by the following 



Proposition 4.6. Let n^,{x.) E A be a local minimizer of the design problem (25) and (m*, w*) 
denote the corresponding minimizing scattering resonance pair. If is nondegenerate and 
|3fJa;^| < p, then 

n+ Q (aw^M^(x)) > 
n_ 3 {auluKx)) < 



where a G C zs defined as in Prop. 4-i- 

Proof. By Prop. 4.3, we may decompose Q = ^2+ U where 

= {x G r2: Ti*(x) = n+} and = {x G ^7: ri^(x) = 

Let A C and 1a be the indicator function for the set A. Local optimality of ra^ implies 
that 



< ^> (^acol j > {aujlul{-K)) > Vx G 



since A is an arbitrary set. The statement for f2_ in (28) follows similarly. □ 



5 Computation of optimal one-dimensional n{x) 

In this section we describe one-dimensional computations that will be used to motivate 
analytical results in section 6. Computationally, we solve a local variant of the optimization 
problem (double-infimum) in (10). Following [HBKW08, KS08], we pick a resonance uj 
and improve the structure n(x), by local gradient methods, to extend the lifetime of that 
particular resonance. Thus we solve the local optimization problem 

Seek local minima of |53a;[n]|, subject to n G A. (29) 

The local minimizers of Eq. (29) which additionally satisfy |3fta)*| < p are local minima of 
Eq. (25). 



5.1 Computational method 

Below, we refer to the forward problem as the computation of the resonances for a given 
n{x) and the optimization problem as the solution of Eq. (29). 
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Figure 2: Left plot: For j = 0, 1, ... ,9, a plot of piecewise constant structure, nl, (blue) 
which minimize (30). Corresponding mode modulus, \ui\ (black) is plotted on the same axes. 
Right plot: Plot of ^ul for j = 0, 1, . . . , 9. 



Forward problem For a one-dimensional {d = 1), piecewise-constant, refractive index 
n{x), the resonances satisfying Eq. (13) are the roots of a nonlinear system of equations 
obtained by imposing transmission conditions at the material discontinuities. In [HBKW08], 
this system of equations is derived and Newton's method is applied for finding the roots 
of this system. This method works extremely well if initialized sufficiently near the desired 
resonance. We initialize Newton's method by either (i) using a finite difference discretization 
of Eq. (13) to form a quadratic eigenvalue problem (QEP) which is solved using Matlab's 
polyeig command [TMOl], or (ii) using the u computed at a previous optimization iteration. 
Other references that derive the transmission conditions at material discontinuities, including 
for two-dimensional, radially symmetric n(x) are [Yeh88, YYM78]. 

Optimization problem. In Prop. 4.1 and App. B, we compute the variation Thus, a 
number of gradient-based methods are available to solve Eq. (29). In [HBKW08, KS08], the 
authors use gradient descent methods. To solve Eq. (29), we have applied a BFGS interior 
point method. 
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Figure 3: The transmission coefficients It-'Xw)! corresponding to the refractive indices ni{x) 
minimizing Eq. (30) for j = 1, . . . , 10. The range of the w-axis is taken to be [0, 23ftt<jf] to 
demonstrate that in the center, \P{^ui)\ ~ 1, while near u = is very small as j 

increases. 



5.2 Computational results 

Let Q = [0, L], L = 1, = 2 and = 1. We define , for j = 0, 1, . . ., to be the scattering 
resonance for which the corresponding mode modulus has j minima^ (see Fig. 2). Using the 
method described in Sec. 5.1, we solve 

nl := argmin | Q^o;-' [n] | . (30) 

for j = 0, . . . , 9. For each j, we plot in Fig. 2 the optimal refractive index ni{x) and the cor- 
responding resonance pairs denoted, (u{,a;f). In Fig. 3, we plot the transmission coefficient 
modulus, |t-'(a;)|, associated with n{{x). The transmission coefficient t{u) is defined by the 
solution of the form 

u(x,uj) = < . for e M 

|t(w)e*^^ x> L 

Remark 5.1. We observe the following: 



^except for mode j = 0, whose modulus has one minima and whose eigenvalue is purely imaginary 
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1. For each j, the optimal n{ is symmetric, i.e., n{ G Asym C A. 

2. For each j, ni(0) = n{{L/2) = ni{L) = n+. 

3. For each j, the optimal ni{x) is roughly periodic with a defect near x = L/2. The 
number of repeated blocks increases as j increases. The width of the repeated structure 
approximately satisfies Bragg's relation as defined in Prop. 6.6. As j increases, for 
even modes, the defect width tends to zero, while for odd modes, the defect width 
tends to approximately twice the Bragg width. These observations will be made more 
precise in Sec. 6. 

4. The values I'^ujH are monotonically decreasing in j, i.e., \'^uji~^^\ < IQ'wf |. 

5. The values are monotonically increasing in j, i.e., the optimal modes u^{x) be- 
comes increasingly oscillatory as j t C)0. 

6 Characterization of locally optimal n{x) in dimension 
one / Bragg relation 

In this section, we fix > 1, n_ = 1, p, and a domain [0, L]. We will use a combination of 
numerical observations from section 5 and analysis to characterize one-dimensional refractive 
indices n^, & A which are 

local minima of r^[n] subject to n E A. (25) 

In this section, we only consider states such that |3?u;^| < p, that is is an interior point 
on the constraint set. Additionally, we take p sufficiently large such that by Lemma 2.3, 
\^u^\ > 0. 

Even / oddness properties of and We begin the discussion with even/oddness 
properties of the optimal refractive indices n^{x) satisfying Eq. (25) and the corresponding 
modes m*(x). The following conjecture is supported by numerical experiments (see Sec. 5.2). 

Conjecture 6.1. There exists G Asym, as defined in (8), such that r^[r;,^] = r^(^). 
That is, the refractive index obtained by minimizing F'' over the set A can be taken to be 
symmetric. 

Motivated by this conjecture, in what follows we restrict consider optimization over the 
constraint set, Asym, of symmetric structures. 

Proposition 6.1. If d = 1 and n G Asym, then any scattering resonance mode satisfying 
the Helmholtz equation (4) is either even or odd with respect to x = L/2. 
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Proof. Let {u{x),uj) be a scattering resonance pair satisfying Eq. (13). Since n{L—x) = n{x), 
the following are also solutions: 

ui{x) = u{x) + u{L — x), U2{x) = u{x) — u{L — x). 

The Wronskian W(ui,M2) = uiU2~u'iU2, is constant in x. The outgoing boundary condition 
implies W(x) = 0. Thus ui and U2 are linearly dependent, which is possible only if either ui 
or U2 vanishes identically. Thus, u{x) is either even or odd. □ 

From the numerical experiments in Sec. 5.2, we observe that modes corresponding 
to locally optimal may be either even or odd (see Fig. 2). It follows that the modes u^, 
corresponding to locally optimal solutions of (25) will be either even or odd depending on 
the choice of p, L, and 

Lemma 6.2. Let n{x) G Agym, i-^- n{x) is symmetric about L/2. Then, if u{x) is a solution 
of the scattering resonance eigenvalue problem (13), then neither u{x) noru'{x) may vanish 
for any x ^ L/2. 

Proof. Let {u{x;oj),co) denote a scattering resonance pair for (13) with n{x) G Asym- By 
Proposition 6.1, we have u{L/2\ w) = or u\L/2] uj) = 0. Suppose there is a ^ G M \ {L/2} 
such that either m(^; u) = or u'(^) = 0. Assume ^ < L/2. If ^ > L/2, the proof is similar. 
Then, u{x] u) is a solution of the eigenvalue problem value problem: 

— dlu{x;uj) = uj'^n'^{x)u{x;uj), ^ < x < L/2 (31) 

with one of the pairs of boundary conditions: 
u{^) = 0, u'{L/2) = or m'(0 = 0, u'{L/2) = or 
u{^) = 0, u{L/2) = or u'{^) = 0, u{L/2) = 



Multiplication of (31) by u{x;uj), integration over the interval [^,L/2], integration by parts 
and using the boundary conditions yields: 

L/2 i>L/2 

\dxu{x;uj)\'^dx = u'^ / n'^{x)\u{x;u)\'^dx (32) 

It follows that cu^ is real and positive. Therefore u is real. However, Theorem 3.1 precludes 
the existence of a real scattering resonance energy. Thus we have a contradiction and Lemma 
6.2 is proved. □ 

The following proposition describes the change in argument of a solution u{x) to the 
Helmholtz equation (4) with arbitrary symmetric index, n{x). 

Proposition 6.3 (Monotonicity of phase). Let {uj,u{-,uj)) denote a scattering resonance 
pair satisfying the Helmholtz equation (4) with n G Asym, i-^-, symmetric about x = L/2. 
Then 

d /"^' 

aiguix) = 2{^u)\Quj\\u{x)\-^ n\z)\u{z)\^dz. (33) 



dx 



L/2 
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// 3^0; > 0, the corresponding resonance state has increasing argument for x > L/2 and 
decreasing argument for x < L/2. The opposite statement holds true for the resonance state 
corresponding to —uJ. In addition, for x G [0,L], 



d 

— arg-ufx 
ax 



< Cu^. (34) 



Proof. A similar proof may be found in [HS86]. By Lemma 6.2, u{x) does not vanish except 
possibly at X = L/2. Thus on the interval [L/2 + e, oo) for small e > 0, we make the Ricatti 
transformation y = ^ and note that '^y = C^ilogu))' = (arg-u)'. We derive a differential 
equation for and show that the solution has the desired properties. Using Eq. (13a) we 
obtain y' = —u'^n^ — y^, which has imaginary part: 

(%)' = 2(3fJa;)|53u;|n2-2(%)(%), 



where we use that = —\'<su\ since scattering resonances lie in the lower half plane. Using 

_d_ 
dx 



the integrating factor exp(log |up) = we obtain the differential equation 



2(3fJa;)|53ci;|r;, |u| . Integrating from + e to arbitrary x we obtain 

argM(x) = |m(x)|-M |M(| + e)|'%(| + e) + 2(SRw)|Sw| n\z)\u{z)\^ dz 



dx 

Taking the limit as e | 0, we find that \u{^ + e)\'^'^y{^ + e) — )■ if u{x) is either even or odd 
giving (33) for x G [L/2, L]. The proof is similar for x G (— oo, L/2]. 

If u{x) is even about x = L/2, Eq. (34) follows from Eq. (33), Lemma 2.5, and Lemma 
6.2 (m doesn't vanish). If u{x) is odd about x = L/2, then |m(L/2)| = and one may use 
L'Hopital's rule to show that 

fX 

-2 / ^2/^M„ /^•m2 



□ 



lim |m(x)|~M n\z)\u{z)\' dz = 0. 



Proposition 6.3 can be used to obtain a significantly more detailed characterization of 
one-dimensional local optima, n^(x), than provided in the general results of Prop. 4.3 and 
Prop. 4.6. In particular, we have 

Theorem 6.4. In spatial dimension one, the optimal refractive index, n^,{x) G Asym, has a 
finite number of transition points between n+ and n_ and is therefore a step function. 



Proof. Our proof is based on the approach taken in [Karll]. Suppose that the optimal 
refractive index, n*(x), has an infinite number transition points, labeled {xj}'jLi, in the 
interval [0,L]. By Proposition 4.6, at each transition point xj, we have '^[a^,ulul{xj)] = 0, 
where is a nonzero constant depending on u^{x) and is in the open lower-half plane. 
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Since the sequence {xj} is bounded, there exists a convergent subsequence, relabeled {xj}'^^, 
which converges to a point X e [0, L/2]. We consider two cases. 

Case 1: Suppose X ^ L/2. By continuity of Uj,(a;), the real sequence ai,uj1u1{xj) con- 
verges to a^wf G M. In addition, the real sequence of difference quotients ai.u'llul^x j) — 
u1{X)\/{xj — X) converges to 2ai,ujlu^{X)u'_^{X) G M. Note that ^ (see Lemma 
6.2). Thus the ratio 

2a,ujlu,{X)u',{X) 
a^u^uliX) 

is a real number, say 2p, implying u'_^{X) = pUi,{X). Thus satisfies a self-adjoint 

boundary-value problem on the interval [X, L/2], contradicting Lemma 6.2. 

Case 2: Suppose X = L/2. Thus a-^ojlul{xj) is a sequence of real numbers converging 
to a-^ulu1{L/2). Equation (33) implies that argM^(x) is monotonic on the interval [0,L/2]. 
Indeed for Sftw* > 0, arg-u+(x) is strictly monotonic. Thus 

I arg [a^ujlul{xj+i)\ - arg [a^ujlul{xj)\ \ > 2n. 

On the other hand, by Eq. (34), we have 

I arg [a^w^M^(L/2)] - arg [a^w^M^(x)] | < Cup\L/2 - x\ 

which is a contradiction. This completes the proof of Theorem 6.4. □ 

Let an optimal refractive index, n^(x), have + 1 discontinuity points, say 1, which 
we denote by {x^j^Q, such that'^ 

X e {xj,Xj+i), j even ^^^^ 
otherwise. 

In Figure 4, we illustrate the relationship between n*(a;) and the sign of Q{aulul{x)) by 
plotting a locally optimal refractive index n^{x) and the quantity Q[aujlul{x)]. 

Since the energy of a mode concentrates where n is large and we expect the energy to be 
concentrated in the center of the interval Q = [0, L], we conjecture the following: 

Conjecture 6.2. ni,{L/2) = n+. 

Conjectures 6.1 and 6.2 imply that the number of intervals with n = n+ is odd. Equivalently, 
the total number of intervals, N = 4M + 1 for M G N. For example, for the refractive index 
in Fig. 4, = 9 and M = 2. We refer to [xo,a;i] as the left-most interval, [xn-i,xn] as the 
right-most interval, [x(^]\f^i)/2,X(^N+i)/2] as the center interval. All other intervals are called 
interior intervals. 

Proposition 6.5. The length of an interior interval Ij = satisfies 

min^e/^ Mx)\'^ 1 2^ 




^Note: we have not ruled out the possibiUty that xo > and xj\j < L. 
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Figure 4: A plot of an optimal n^(x) — 1 in black and '^[auj1u1{xy\ in red. Note that Xq = 
and ^{aulul^Q)) 7^ 0, but the optimality condition (28) is satisfied. This optimal structure 
and mode were obtained by solving Eq. (30) for j = 4. 



Proof. Corollary 4.6 implies that the argument of m*(x) changes by exactly 7r/2 on an interior 
interval Ij = [xj,Xj+i]. Using Eq. (33), we compute 

— = 2|9fta;^| / |u^(x)|^^ / nl{z)\Ui,{z)\'^ dz dx 

2 Jxj Jl/2 

< 2\^uJ l^LoJ I \Ui,(x)\~'^ / nl(z)\uJz)\'^ dz dx 



L/2 



\^uj^\ [^^(O)^ / |M^(x)|~^dx < \^u^\ |m^(0)P \xj+i — Xj\ max |-u^(x)|~^ 



where we used Eq. (15) and the symmetry of n{x) and |m^(x)|. Now by Lemma 6.2, [^^(a;)! > 
for all X & Ij, from which Eq. (36) follows. □ 

6.1 Locally optimal refractive indices and the Bragg relation 

Denote a{x) := '^{auluKx)). From Proposition 4.6, 

cr(xo) = argM^(xi) - argM^(xo) = -. (37) 

We remark that numerical experiments suggest that cr{xQ) ~ (see Fig. 4) and that cr(a;o) — )■ 
as n+ or L — )■ CX3. For the idealized limiting case we have: 

Proposition 6.6 (Bragg relation). Assume a locally optimal mode satisfies (37). Then, the 
intervals on which n^(x) is constant have the following properties: 

1. The length of all non-center intervals, i.e., Ij = for all j 7^ (A^ — l)/2 is 

given by d± = \ n±^^uj^,\ ~ \^eff,±, where Xeff,± is the "effective" wavelength of u^, in the 
medium n±. 
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2. Thus, a full "period" of the repeated structure is S = + = 2n^'\^5~\' "^here 
Uh = '^{fH^ + ^2 ^)^^ is the harmonic mean of and 

3. Additionally, the length of the center interval is less than 2d+. 

Since the numerical evidence is that (37) is satisfied only asymptotically for large Uj^ or L, 
this suggests the following: 

Conjecture 6.3. The widths of the non-center intervals, 
w r^ 1 27r 1 

d±{n+,n^, L) -— — , — = -Aefr± as ri+ — )■ oo or L — )• oo. 

Proof of Proposition 6.6. Write Ui, = ujr + iuj. The solution on any interval Ij = [xj, x^+i] 
where n^(x) = n± can be written 

Ui,{x) = ttj cos[a;*n-|-(a; — xj)] — if3j/n± sin[a;*?T,-|-(x — xj)] (38) 

where aj and l3j are chosen such that ) and Pj = {i/uj)u'^{xj). We record the 

following formulas for use in the argument below: 



u 



■l(^) = —ajOJ-^n± sin[c<;^n-|-(x — Xj)] — iijj^,f3j cos[u^,n±{x — Xj)] 



sm = cosh , cos = —t smh 



\2ujrJ \2ujrJ ' \2ujrJ \2ujr 

2ur 



u[{xj + d±) = —ajUJi,n± cosh [ - — - ) — u^Pj sinh 



(39) 



Without loss of generality, let > and assume that on the first interval, Iq = [xq, xi], 
tto = /3o = 1 which satisfy the outgoing boundary conditions for x < xq. We now show that 
as X increases, there are two alternating cases which define n^{x) on the interior intervals. 

Case 1 Suppose that on the interval Ij = [xj, xj+i], aj, Pj G M and = This is the case 
for the first interval, Iq. On an interior interval, by Prop. 4.6, the point xj+i is that point at 
which argM^(x) has increased by 7r/2. On the first interval, this is precisely assumption (37). 
Since ^Ui,{xj) = 0, an increase of argM^(a;) by 7i/2 on Ij is equivalent to 3ftuj,(a;j+i) = 0. 
The value of x such that 9fJM*(a;) given in (38) next vanishes is Xj+i = Xj + d+. Using (39), 
u'{xj+i) = ujj where 7 G M. Thus, the coefficients a^+i and on the next interval, Ij+i, 
are purely imaginary. 



Case 2 Suppose that on an interior interval, Ij = [xj,Xj^i], n^{x) = n_ and aj and Pj 
are purely imaginary, i.e., aj,Pj G iM.. By Prop. 4.6, the point Xj+i is defined to be the 
first value for which '^u^,{xj+i) = 0. Using, (38), we find that Xj+i = Xj + d-. By (39), 
u'{xj+i) = iu'y where 7 G M. Thus the coefficients Oj+i and Pj+i on the next interval Jj+i 
are purely real, which is precisely Case 1. 
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Thus, the intervals alternate between Cases 1 and 2 until the center interval, I[n~i)/2 = 
[a;(Ar_i)/2, X(Ar+i)/2]. is reached. Here, by Prop. 6.3, argu*(x) increases on [x[n~i)/2iL/2] and 
decreases on [L/2^X(n+i)/2\- By Prop. 4.6, argM^(x) changes by less than 7r/2 on each of 
these half intervals, which implies that the width of each of these half intervals is less than 
d+. □ 



7 Discussion 



In this article we have investigated the question of finding compactly supported structures, 
defined in terms of a refractive index, n(x), for which the Helmholtz equation has very 
long lived scattering resonances in dimension d =1,2,3. Here n varies in an admissible set 
of refractive indices ^ (fi, n_, n+, ), defined in terms of the support of the structure 
and upper (n+) and lower (n_) pointwise bounds for n. We now summarize, based on 
our analytic, asymptotic, and numerical studies, further expectations, beyond what we've 
established, for the optimal one-dimensional structure and its corresponding resonance of 
minimal imaginary part, l^^w*! = F^, 

1. We expect that the optimal resonance such that r^(^) = I^Jw^l is achieved for 
a resonance with real part approximately imposed upper bound: 3fJa;* ~ p. We also 
expect its imaginary part to be exponentially small in p: '^u^, ~ 74exp(— where 
A and B are constants dependent on A. 

2. Furthermore, we expect that the number of times the optimal refractive index, n*(a;), 
alternates between n+ and ri- on [0,L/2] is approximately M = (N — l)/4 ^ ^ = 

3. Simple models [HW08] suggest that \^oo^,\ ~ n^^'' and thus F^ ~ n^^^^^^'^^ . 



Finally, we remark on a number of interesting future questions and directions. In the 
present work we have established that classical one-dimensional Bragg structures arise as 
optimal one-dimensional compactly supported dielectric cavities with finite material contrast 
in the limit that the support or material contrast tends to infinity. Future work will include 
study of the deviation of locally optimal from the Bragg relation defined in (6.6) for 

finite support and material contrast. 

Our work establishes the existence of a solution to the optimal design problem in two 
and three spatial dimensions {d = 2,3). It would be interesting to study the limit of large 
support or large material contrast for c? = 2 or 3. Perhaps this leads to a higher-dimensional 
variant of Bragg structures. 

We conclude mentioning that the measure of resonance lifetime to be optimized is often 
chosen to suit the particulars of the problem. We have chosen to maximize: r = |53cj|~^. 
Other choices include the quality factor, Q = which measures loss per unit oscillation 
cycle, and the Purcell factor, oc Q/V, where V denotes the "mode volume". The Purcell 
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factor is particularly important in applications where a strong light-matter interaction is 
desired. Maximization of such quantities would be a natural and interesting extension of 
this work. 



A Examples: resonances for radially symmetric cavi- 
ties in M^, ti= 1,2,3 

In this subsection, we discuss the resonances for a simple example in dimension d = 1, 2, 
and 3. We take f2 = {x: |x| < a} and n(x) defined by: 



n(x) = <^ " ' ' (40) 




where hq > 1 and a > are constants. We now consider the scattering resonance problem 
(5); see also [Noc97, Ch. 8] and [DBD+08]. 

A.l Dimension d=l 

Imposing outgoing radiation conditions we find that 

{Ae-'^"" X < -a 

De''^'' x> a . 

Imposing continuity of u and d^u at x = ±a yields a 4 x 4 linear system of equations 
for the constants A, B, C, and D, whose non-trivial solvability requires the vanishing of 
a determinant. Resonances are the values w G C for which this determinant vanishes: 
^-2iuj{no-i)a |-g4iajnoa^^^ _ _ _|_ -|^^2j ^ whosc solutious are giveu by 



nm 1 

CUm = i log 



no + 1 



riQ-l 

The resonances in Eq. (42) for no = 2, a = 1 are plotted in Fig. 5(left). 
A. 2 Dimension d=2 

Outgoing solutions bounded at the origin are given by 



m G N. (42) 



u{r, 6) 



AJi(nQUjr)e^^^ r < a 
BHf\ujr)e'^'^ r > a, £ G Z 



Imposing continuity of u and drU at r = a leads to a system of two linear homogeneous 
equations for A and B. This equation is solvable if and only if 

j£{noUja)Hl^^' (uia) — nQHf\uja)J'^{nQUJa) = 0. (43) 
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Figure 5: Helmholtz resonances for index of refraction given in Eq. (40) in 1-, 2-, and 3- 
dimensions computed using Eqs. (42), (43), and (45). In each dimension, tiq = 2 and a = 1. 



We numerically solve Eq. (43) for no = 2, a = 1, and £ = 0, . . . , 9 and plot the resonances 
in Fig. 5 (center). For each angular momentum, £ > 0, there is an infinite sequence of 
resonances, {w^j}j>i, in the lower half plane. Fixing I and letting j tend to infinity, these 
resonances are seen to have the asymptotics: 



vr 
riQa 



j + i + 



— 1-. 



2noa 



log 



no + 1 



no 



J » 1, 



(44) 



as derived in [Noc97] . For fixed j and angular momentum, i, tending to infinity the sequence 
of resonances approaches the real axis exponentially [Noc97, DBD+08]. Note that Uqj has 
multiplicity 1 and uej has multiplicity 2 for j > 1. The mode plotted in Fig. 1 (center) 
corresponds to the 2-dimensional resonance with £ = 6 and j = 1. 



A. 3 Dimension d=3 

The analysis is analogous to the case of spatial dimension 2. Solutions bounded at the origin 
and outgoing at infinity are given by 



u{r, 6, < 



Aji{noUjr)Yl^{9,(l)) r<a 
Bhf\ur)Y^{e, 0) r > a, |m| < £ G N . 

Imposing continuity of u and drU at r = a leads to the solvability condition 



jf{noUja)hf^' (uia) — nQh\^\uja)j[{nQUja 



0. 



(45) 



We numerically solve Eq. (45) for no = 2, a = 1, and £ = 0, . . . , 9 and plot the resonances 
in Fig. 5(right). For fixed angular momentum, £, and j tending to infinity, the scattering 
resonances follow the asymptotics 



TT 



no a 



] + 



2noa 



log 



no + 1 



no - 1 



m >> 1. 



(46) 



For fixed j and large angular momentum, i, an asymptotic analysis involving spherical Hankel 
functions, analogous to that in [Noc97, DBD"'^08], shows that the uej tend exponentially 
toward the real axis. Note that uej has multiplicity 2i + 1. 
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B Calculation of the variation of a scattering reso- 
nance, cj, with respect to the refractive index, n(x); 
proof of Proposition 4.1 

The scattering resonance problem (5) can be written 

m(x;w) = uj^R^[u]{x;uj), where R^[u]{^;uj) := / G'(x,y,w) m(y) u{y;uj)dy (47) 

Jn 

and m(x) := n'^{x.) — 1 is non-negative with support in Q. 

We consider a perturbation of G ^ of the from ?t,(x) i— )■ n(x) + 5n(x) where 6n{x.) G 
L°°(]R'^) with support in Q. Denote variations with respect to n using a prime and the 
variation of G{'x, y, u) with respect to u by (^(x, y, w). For x G we take variations of (47) 
to obtain 

(id — uj'^Rui) u{'k; u) = 2ujuj' Ruj (x; uj) + uJ^ I (^(x, y , uj) J m(y ) u(y ; w) dy 

Jo. 

+ G(x,y,u;) m'(y) M(y;u;)dy. 

Jn 

Denoting {f,g) = fg and taking the inner product with a general v G L'^{Q), we obtain 
{v, {ld-uj'^R^)u') =2coco'{v,R^[u]{-;co)) (48) 




+ UJ u / t;(x) G'(x,y,a;) m(y) M(y;w)dydx 



n Jn 



+ / t;(x) G'(x,y,a;) m'(y) M(y;a;) dydx. 

Jn Jn 

Using the adjoint operator 

KlvK^;^^) = m{x) / G{x,y,uj)v{y)dy, 
Jn 

the left hand side of (48) can be rewritten 

{v, {id - Lo'^Ru) u) = / I;(x)m(x) dx - / / t;(x)G'(x, y, w) m(y) ^'(y; w) dy dx 

Jn Jn Jn 

= {{ld-uj'R:)v,u'). 



Setting v{-k) = v{-k;u) = u m{x.)u{-K;uj), we compute 

(Id - uj'^Rl) t;(x; u) = uj'^m{x) |^u(x; u)-uj'^ J ^(x, y, u;)m(y)M(y; u) dy 



u^m{x){ld-uj^R^)u{x;uj) =0. (49) 
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Therefore, for the above choice of ^(x;^), the left hand side of (48) vanishes. The first and 
third terms on the right hand side of (48) are evaluated to be: 



{v,R^[u]) = (u'^mu^uj ^m) = (m u^u) 



and 



(50) 



(51) 



/ / t;(x) (^(x, y, uj) m'(y) u{y; u) dy dx = (m, m'u). 
Jn Jn 

Therefore (48) reduces to 



u' i 2uj{u, m u) + u'^ / m(x) m(x; tu) (^(x, y , u) m(y ) ^(y ; u) dy dx 1 = — (m, m' u) 



5(jj 



n Jn 



or, equivalently, recalling that m' = (n^)' = 2n6n, 



UJ 



where 



a = 2u{u{-;u),m u{-;u)) + u / / m(x) m(x; w) ^(x, y, w) m(y) n(y; w) dy dx. 

Jn JQ 



Thus, 



with 



u[n + 6n] — u[n] — ( —,5n 

on 



6n 



—2 a UJ nu{-;uj) , 



proving (26). 

In one dimension, the explicit formula for the Green's function given in (3) may be used to 
obtain the first equality in (27). This formula agrees with [HBKW08, Eq. (19)]. The second 
equality in (27) can be obtained using the identity cu^ / n'^u^ = J ul — za;[u^(0) + u^{L)], 
obtained by multiplying (13a) by u and integrating by parts. 



C Bragg's relation maximizes the spectral gap to 
midgap ratio 

Consider a one-dimensional n{x), infinite in extent, which is periodic with period d, i.e., 
n{x + d) = n{x), and has alternating layers, i.e., 

\ rii < X < b 
n{x) = < 

\n2 b < X < d. 
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Figure 6: (left) The dispersion relation (52) for a quarter wave stack where ^ = • There 
is a spectral band gap centered at w = ^ of width ^ sin~^ (right) For ni = 1 and 

n2 = 2, we let 6 = and plot 7 vs. -R(7), the spectral gap to midgap ratio as defined 

in (55). We find that R{'^) is maximal precisely when 7 = 1, the structure defined by the 
Bragg relation. 



It is shown in [Yeh88] that the wave equation (1) has a solution of the form v{x,t) = 
e^^^^~^^^Uper{x) where Uper{x + d) = Uperix) if uj and k satisfy the dispersion relation 

cos{kd) = cos{unib) cosun2{d — b) — - ( — -\ — -\ sin(a;ni6) smun2{d — b). (52) 

2 \ni n2 J 

The Bragg relation is defined 

h M M ^ ^ 1 2n rihd rihd 

mb = n2[d - b) = - — =^ d=- , b=—-, [d-b) = —- (53) 

4 oj Irih OJ rti 2 n2 2 

where Uh = 2(n|f^ + n2^)^^ is the harmonic mean of rii and n2. The width of each layer is 
a quarter wavelength, i.e., the phase of the wave changes by 7r/2 in each layer. For a Bragg 
structure, the dispersion relation (52) simplifies: 

co.(fcci) = l-l^^^i±^sin^ (^). 
^ ' 2 nin2 \ 2 ) 

Thus, the dispersion relation for a Bragg structure has a spectral band gap centered at 
u = of width -1^ sin~^ \ni-n2\ ^ rpj^^ dispersion relation with n^ = 1 and no = 2 is plotted 

anil dnii ni+n2 ^ j- ^ r- 

in Fig. 6 (left). 

The Bragg relation given in Eq. (53) can be interpreted to state that constructive in- 
terference occurs when the path length of a reflected wave is equal to (a multiple of) the 
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wavelength. In fact, the Bragg structure is optimal in the following sense. For fixed ni, n2 
and d, let 

bil) = l]— = l]x,s, 7e(0,l) (54) 
4 Lorii 4 

define a class of periodic layered structures. The Bragg structure is defined by 6(1). Let Ui 
and U2 be the left and right edges of the first spectral bandgap. For a given layered media, 
define the spectral gap to midgap ratio by 

R := (55) 

In Fig. 6(right), we plot the the spectral gap to midgap ratio for the class of devices defined 
as in (54). The maximal gap to midgap ratio occurs for the Bragg structure (7 = 1). Lastly, 
we note that the Bragg structure does not maximize the spectral gap width, given by |ci;2— cJil 
[Ostll]. 
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